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Abstract 

We investigate the linear and nonlinear optical responses of dilute anisotropic 

networks using the Green's-function formalism (GFF)[Gu Y et al. 1999 Phys. 

Rev. B 59 12847]. For the different applied fields, numerical calculations 

indicate that a large third order nonlinear enhancement and a broad infrared 

absorption arise from the geometric anisotropy. It is also shown the overlap 

and separation between the absorption peak and nonlinear enhancement peak 

when the applied field is parallel, perpendicular to the anisotropy respectively. 

In terms of the inverse participation ratios (IPR) with q = 2 and spectral 

distribution of optical responses, the results can be understood. 
PACC: 7430G, 8160H 



1 



Typeset using REVT^ 



I. INTRODUCTION 



Recently, optical properties of composites have attracted great interest. In order to open 
new possibilities in the information processing and transmission, a large optical nonlinearity 
may be desirable 0. Composite materials, especially, small metal particles embedded in a 
dielectric host and metal clusters on the nanometer scale, exhibit a strong nonlinear optical 
enhancement through the inhomogeneous local-field and geometric-response effect frl. It 
is also known that these composites give rise to an anomalously large absorption in the 
infrared spectrum . To analyze the optical responses, it is more convenient to adopt the 
spectral representation ^. For various anisotropic composites, the spectral density was cal- 
culated as a function of volume fractions p\\ and p± by the effective-medium approximation 
(EMA) @§§. It was found a large nonlinear optical enhancement, as well as the separation 
of the absorption peak from the nonlinear enhancement peak. By EMA, the fluctuation of 
local field has been averaged out. So it is difficult to find the physical origins of the opti- 
cal enhancement and the separation of optical peaks. In this connection, Green's-function 
formalism (GFF) was developed to deal with the optical responses of the arbitrary-shaped 
metallic clusters embedded in the infinite dielectric networks at the quasistatic limit 00. 
By the formalism, the resonance spectrum and local field distribution for each eigenmode 
can be analytically obtained. The aim of this paper is to investigate the linear and non- 
linear optical responses throuth the dilute anisotropic networks in view of the local field 
distribution. 

In the following, a binary dilute anisotropy network is considered. We use the random 
generator to produce the geometric anisotropy. The metallic bonds parallel to the applied 
field are assigned with the fraction p\\, and metallic bonds perpendicular to the applied 
field with p±. In the dilute systems, the condition, p\\ x p^ = 0.01, is satisfied. When we 
change p\\ from 0.1 to 0.8, the corresponding p± varies from 0.1 to 0.0125. In Section II, 
the Green's functions for the effective linear response eg and effective nonlinear response Xe 
are derived. In Section III, the inverse participation ratios (IPR) with q = 2 are used to 
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represent the localized and extended eigenstates in anisotropic networks. In Section VI, the 
spectra of the absorption and the third order nonlinear enhancement are illustrated. When 
the applied field is parallel to the anisotropy of the networks, the overlap of the absorption 
peak with nonlinear enhancement peak is found within the interval 0.1 < p\\ < 0.8. In 
contrast, for the perpendicular applied field, the separation of the absorption peak from 
nonlinear enhancement peak is enhanced when the anisotropy p\\ is increased. These results 
are explained in Section III, numerically confirmed in Section VI, and concluded in Section 
V. 

II. EFFECTIVE LINEAR AND NONLINEAR RESPONSES 

Consider an infinite binary network as shown in Fig. 1, where the impurity bonds with 
admittance ei are employed to replace the bonds in an otherwise homogeneous network of 
identical admittance €2- The admittance of each bond is generally complex and frequency- 
dependent. All the impurity bonds construct the clusters subspace. When resonance hap- 
pens, the potential of the jointing points can be computed by the GFF. Instead of the point 
source 6^^, the source term is replaced by px at the point x. Hence in the subspace, V is 
a linear combination of right eigenvector R^s of Green's-matrix M, 

V=Y: TT E iLn,y E P^'Gy,.')Rn. (l) 

n=l ^2 1,* i>n) ygc x' 

And for the site x(a;i,a;2) outside the cluster, becomes 

K = E Px'Gx,x' + E V E (^"-y E P^'^y^^') E ^x,z^2n,z. (2) 

x' n=l ^2(,S Sn) ygc x' zeC 
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FIGURES 



FIG. 1. Schematic diagram of a cluster (shown in thick lines) embedded in an infinite network 
In the uniform field Eq along 1 direction, V reads 

V=t ,(f\A J:Ln,yyi)Rn. (3) 
n=l ^21-5 - Sn) 

For a binary network with N x N square lattices, in the quasistatic limit, the displacement 
D of the ith. bond is related to the local field E by the relation D = eiE + Xi|E|^E 
where is the dielectric constant and Xi is the third order nonlinear susceptibility of the 
bond. ei,Xi are set to be ei, xi for the impurity bonds, and €2, X2 for the matrix bonds 
respectively. By the finite difference transformation, the effective linear response along the 
applied field is 

(x,y)eC " (x',y')^C " 

with the magnitude of field Eq = 1 and lattice constant a = 1, where (x, y) G C means 
that X and y are the nearest neighbors in the same cluster while (x', y') ^ C means x' and 
y' are the nearest neighbors but not in the same cluster. When we employ A^(A^ + l)e2 = 
Sx,y^2(K; — V^), 62 — cau be cxprcssed by the points within the clusters subspace as, 

iV(iV + l)e2 - ee = Yl (e2-ei)(K-^y). (5) 
(x,y)ec 

Here ei/e2 = (s — so we have 

e 1 ^ ^ ^ 

N{N+1)-^ = Y.- r(EWl) E iXl-yi)iRn,.~Rn,y). (6) 

£2 n=l ^2[S - Sn) (x,y)gC 

Because x and y are neighbors, the values of Xi — yi can be only +1 or —1. Physically, the 
absorption along the applied field is larger than that along the perpendicular direction of 
the applied field. In the following calculations, ee always represents the absorption along the 
applied field. For simphcity, we let Lf, = Eyec Ln,yyi , Rh = Ex,yGc(a;i - l/i)(-Rn,x - -Rn,y), 
and 7„ = L^R^, the imaginary part of corresponding to the absorption is 



-Im{e,) = Imt-j^y (7) 

n=l ^2l'5 - Sn) 

7„, referred to as the cross section, obeys the sum rule 0, 

E^n = N, (8) 
n 

where Nh is the number of horizontal bonds along the applied field. 

Since the local fields are determined completely, our formulas can be used to calculate 
the nonlinear response. By relating the total electrostatic energy to the effective coefficients, 
the third order nonlinearity x'e is defined by 

D(x)E(x) = V^KE2 + xir], (9) 

where E = (l/V) /„E(x)d^x is the space averaged electric field. For the infinite networks, 
E = Eq = 1 and the nonlinear response function is given by S. 

X'e = ^J2x^\SViW, (10) 

where the summation is over all bonds and 6vi is the(general complex) potential difference 
across the bond i and for two dimensional case, V = t^. When all the bonds have nonlinear 
term, i.e., Xi = X2 = 1-0, with the lattice constant a = 1, the effective nonlinear response is 
written as 

Xe = V^xl = E l^x - Vy\\V^ - V^f. (11) 
(x,y) 

The third order nonlinearity Xe is expressed as the summation of the forth moment of local 
electric field. So the fluctuation of local electric field enhances the nonlinear optical responses 
well. 



III. INVERSE PARTICIPATION RATIOS WITH Q = 2 



First, the IPR of eigenvectors of Green's-matrix M in the eigensystem with Ug jointing 
points are defined. The nth normalized right eigenvector is 
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Rn — {Rn,li ■ ■ ■ 1 Rn,ii ■ ■ ■ Rn,ns\ (-^^) 

with < i?„ >= and < Rn >= 1. The IPR of Rn is written as 111, 

lPR{Rn)=J2Rn/'- (13) 

i=l 

Here q = 2. The calculation of IPR(-R„) is limited in the nontrivial eigenstates. The number 
of IPR(-R„) is equal to or less than n^. In the above Eqs. (1) and (2), it is found that the 
right eigenvectors of M are closely related to the local fields of the impurity cluster in the 
subspace. So the IPR can be used to represent the localization of the eigenstates. The IPR 
amplify the profiles of eigenstates, namely, the localized states become more pronounced and 
extended states become smoother. Hence, the larger values of IPR are always corresponding 
to the stronger optical responses. This will be verified in the following sections. 

To consider the size effect, the IPR in the dilute isotropic composites with p\\ =0.1 and 
Pj_ = 0.1 are shown in Fig. 2. In this figure, two samples are in size 30 x 30 and 60 x 60. 
The peaks represent the localized states and the valleys correspond to the extended states. 
Comparing the distribution of peaks and valleys in two samples, we find that the size effect 
is not very obvious. It is also seen that the density of states are larger around s = 0.5 than 
that around s = 0.0 or s = 1.0. The localized states incline to accumulate s = 0.0 or s = 1.0, 
while, the extended states are near s = 0.5. So there exists a duality about s = 0.5 0. The 
high values of IPR imply the strong optical responses at s = 0.0 or s = 1.0. 



FIG. 2. Size effect of IPR of right eigenvectors for two samples: 30 x 30 and 60 x 60. Here 
P\\ =Pi = 0. I, pj_=p2 = 0. 1. 

Fig. 3 displays the IPR of the dilute anisotropic systems. We have known that the 
resonance spectrum is very sensitive to the microstructure ©. For each case, comparing 
two samples, we see that the peaks of IPR are very stable though their microstructure 
is completely different. So in the following, for the definite parameters p\\ and p±, it is 
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reasonable to investigate the optical responses by only one sample. When the anisotropy 
is increased from Fig. 3(a) to Fig. 3(g), the localized states incline towards s = 0. or 
s = 1. 0. The distributions of localized and extended states have roughly the same feature 
as that in the isotropic case, i.e., we find more localized states at s = 0.0 and s = 1.0 than at 
s = 0.5, but denser at s = 0.5. The distribution of IPR is caused purely by the morphology 
of the sample. In the next section, we will discuss how the different applied fields act on the 
optical responses when the geometric anisotropy is increased. 

FIG. 3. IPR of right eigenvectors of the dilute anisotropic systems. Sample 1 and sample 
2 are randomly chosen. In (a) and (b), = pi = 0. 2, p± = P2 = 0. 05. In (c) and (d), 
P\\ = Pi = 0. 4, p± = P2 = 0. 025. In (e) and (f), = pi = 0. 5, p± = P2 = 0. 02. In (g) and (h), 
P||=pi = 0. 8, p^=p2 = 0. 0125. 

IV. SPECTRA OF ABSORPTION AND THIRD ORDER NONLINEAR 

ENHANCEMENT 

The optical responses can be properly described by the absorption spectrum 7„ and the 
third order nonlinear enhancement spectrum as shown in Fig. 4 and Fig. 5. When the 
applied field is parallel to the direction of anisotropy, it affects the optical properties well. 
The envelopes of the absorption and third order nonlinear enhancements become narrower 
and the intensity is more enhanced with the increasing anisotropy. The main optical peaks 
are ranged in the interval s G [0.5,0.6], which correspond to the dipolar approximation. 
In this figure, we can not find the separation of the absorption peak from the third order 
nonlinear enhancement peak. The results of single sample conflict with those of previous 
works '0. The reason is that the EMA has averaged out the fluctuation of local field, 
which plays the central role for the effective optical responses. At s = 0.5, the density of 
localized states are much smaller than that of extended states. The optical properties of 
the specific sample are determined by its own eigenstates, not by "effective medium" . So it 
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is not reasonable for all of the cases to express the optical properties of one specific sample 
only by the same macroscopic parameters p|| and p±. In Fig. 5, when the applied field 
is perpendicular to the anisotropy, the optical responses of the dilute anisotropic samples 
are plotted. It is seen that the peaks of absorption are red-shifted while the peaks of the 
third order nonlinear enhancement are blue-shifted when the anisotropy is added. There is 
a larger peak separation with increasing geometric anisotropy. For all the applied fields, it 
is natural to find a large third order nonlinear enhancement and a broad infared absorption 
because the localized states exist in the whole resonant area. Only for the perpendicular 
applied field, the separation of optical peaks is found. 



FIG. 4. Spectra of the absorption and the third order nonlinear enhancement of the dilute 
anisotropic systems for the parallel applied field. In (a) and (b), py = pi = 0. 2, p± = P2 = 0. 05. 
In (c) and (d), = pi = 0. 4, p± = P2 = 0. 025. In (e) and (f), py = pi = 0. 5, p± = P2 = 0. 02. 
In (g) and (h), py = pi = 0. 8, p^^ = p2 = 0. 0125. 



FIG. 5. Spectra of the absorption and the third order nonlinear enhancement 

of the dilute anisotropic systems for the perpendicular applied field. In (a) and (b), 
P\\ = Pi = 0. 2, p± = P2 = 0. 05. In (c) and (d), p\\ = pi = 0. 4, p± = P2 = 0. 025. In (e) 
and (f), py = pi = 0. 5, p± = P2 = 0. 02. In (g) and (h), py = p^ = 0. 8, p± = P2 = 0. 0125. 

Then, the Drude free electronic model is employed to calculate the optical properties of 
the dilute anisotropic composites. The admittance of the impurity metallic bonds is 

. (14) 

where uj is the plasma frequency, and 7 a damping constant. For metal, the plasma frequency 
Up 10^^, being in the ultraviolet. We choose 7 = O.la;^, which is the typical value of a metal, 
and €2 = 1.77, which is the dielectric constant of water for model calculations. The range 
of optical responses is a;/a;p e (0, 1). 



Optical responses of various anisotropic networks are illustrated in Fig. 6. The nonlinear 
enhancement x^^^(or Xe) is divided by 100 for the comparison with the linear absorption. 
The direction of anisotropy in Figs. 6(a), (c), (e) and (g) is parallel to the applied field. 
The peak at u/ujp = 0. 6 corresponds to the dipolar approximation of the whole system. 
In these figures, the separation of absorption peak and nonlinear enhancement peak is not 
found, as observed in the above Fig. 4. Note that this result is different from that of the 
previous works ffl0. However, when the applied field is perpendicular to the anisotropy, 
as shown in Fig. 6(b), (d), (f) and (h), in each of which, both the separation of linear and 
nonlinear peaks, and the third order nonlinear enhancement are seen. 



FIG. 6. Linear and nonlinear optical responses of the dilute anisotropic systems by Drude 
model. Para. field means that the applied field is parallel to the anisotropy while perp. field is that 
the applied field is perpendicular to the anisotropy. Here pi is corresponding to in the text and 
P2 is p±- 

V. DISCUSSION AND CONCLUSIONS 

In this paper, optical responses of dilute binary anisotropic networks are investigated 
numerically by the GFF. Green's functions for the effective linear and nonlinear responses 
are derived. The IPR with q = 2 are given to describe the localization and extension of 
eigenstates. The spectra of the absorption and the third order nonlinear enhancement are 
exhibited for the different applied fields. We also compute the optical responses for the 
Drude model. For the dilute anisotropic networks, we conclude that: 

1. The size effect of optical responses is not obvious. 

2. The IPR with q = 2 can be used to represent the localization and extension of eigen- 
states. 
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3. The wide optical absorption and large nonlinear enhancement are caused by the geo- 
metric anisotropy. 

4. The peaks of the absorption and nonlinear enhancement overlap when the applied field 
is parallel to the anisotropy. While, for the perpendicular applied field, the absorption 
peak is separated from the nonlinear enhancement peak. 

5. The structure sensitivity does not affect the main optical properties, but the details 
of optical responses. 

/^From above 2, 3, 4 and the previous works ffl0, geometric anisotropy is the main 
reason to enhance the optical responses. For a specific sample, when the different fields are 
applied, the positions of optical responses do not vary, but the local field distributions are 
quite different. Therefore the optical properties of resonant composites are determined by 
the geometric anisotropy, as well as by the applied sources. 
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